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Abstract 

m 

Explicit solutions of the classical Calogero (rational with/without harmonic confin- 
ing potential) and Sutherland (trigonometric potential) systems is obtained by diago- 
nalisation of certain matrices of simple time evolution. The method works for Calogero 
& Sutherland systems based on any root system. It generalises the well-known results 
by Olshanetsky and Perelomov for the A type root systems. Explicit solutions of the 
(rational and trigonometric) higher Hamiltonian flows of the integrable hierarchy can 
be readily obtained in a similar way for those based on the classical root systems. 



1 Introduction 

The classical and quantum integrability/solvability of Calogero- Moser systems pP [3j man- 
ifests itself in many guises; the existence of Lax pairs and/or Dunkl operators, algebraic 
linearisation, quadratic algebras, associated integrable spin chains, 'quantised' classical spec- 
tra, etc. Among them, a most intuitive understanding of solvability/integr ability is provided 
by the fact that explicit solutions of the classical equations of motion are obtained by diag- 
onalisation of certain matrices having trivial time evolution as shown by Olshanetsky and 
Perelomov [Hj for the rational and trigonometric potential cases [7]. Their results are for the 
systems based on the A type roots. Here we will show that the same results hold universally 
for systems based on any root system. To be more precise, for the rational potential (with- 
out/with harmonic confining potential) cases, the diagonalisation method works for any root 



system including the non-crystallographic ones. For the trigonometric and hyperbolic poten- 
tial cases, we show that the explicit diagonalisation method holds for any crystallographic 
root system based on the universal Lax pair jH] . A simpler form of explicit diagonalisation 
is provided by the minimal Lax pair j^J E], which exists only for those based on the A, 
D, Eq and E 7 root systems. The basic idea of the explicit solution method is very closely 
related to the notion of algebraic linearisation proved universally for any root system by 
Caseiro-Frangiose-Sasaki jTUj. We will follow the notation of JU] throughout this paper and 
eq.(a.b) of this paper will be cites as (Ia.b). Explicit solutions in terms of diagonalisation 
is readily obtained for the higher (rational and trigonometric) Hamiltonian flows belonging 
to the integrable hierarchy. This works, however, only for those based on the classical root 
systems, A, B, C and D. The conventional Lax pair in terms of the set of vector weights 
(A, C and D) or the set of short roots (B) is indispensable. 

This paper is organised as follows. In section two, the historical background and the 
logical structure of the Calogero-Moser systems necessary for the present paper are briefly 
reviewed. The Hamiltonian and the universal Lax pair with rational potential are introduced. 
The explicit integration in terms of diagonalisation is achieved by relating the Lax pair ma- 
trices L and M to a matrix W of the same size with trivial time evolution, W = 0. Section 
three is devoted to the explicit solution of the systems with rational plus the harmonic con- 
fining potential. In section four, we show the explicit solutions of the Sutherland systems, 
which have trigonometric/hyperbolic potentials. In §4.1 a simple form of explicit diago- 
nalisation is obtained by re-interpreting the formulas of the algebraic linearisation method 
developed by Caseiro-Frangoise-Sasaki [TQ] . This is based on the minimal Lax pair, which 
exists for those based on A, D, E§ and E 7 root systems. A general treatment of explicit 
integration of the Sutherland systems in terms of the universal Lax pair is provided in §4.2. 
This applies to any crystallographic root system. Sections five and six are devoted to the 
problem of explicit integration of the higher Hamiltonian flows of the integrable hierarchy. 
The rational potential case is discussed in section five and the trigonometric case in section 
six. The final section is for summary and comments. 

2 Rational Potential 

The integrability of the Calogero-Moser systems has a long history. Firstly, various types of 
the integrable potentials are recognised; starting from the Calogero model pQ with rational 
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(1/g 2 ) plus a harmonic confining (q 2 ) potential followed by the Sutherland model with 
a trigonometric (1/ sin 2 q) potential. Then the pure rational potential (1/q 2 ) |S] and the 
hyperbolic (1/ sinh 2 q) jl] and the elliptic (p(q)) potentials are added to the list of the 
integrable potentials. As seen in many other subjects in mathematical physics, the quan- 
tum groups, the integrable spin chains, Yang-Baxter equations, etc., the elliptic case is the 
generic one giving all the rest in various degenerations. However, each degenerate case, the 
rational, trigonometric and the hyperbolic, has its own special properties and merits not 
shared by the more generic ones; for example, the algebraic linearisation ^U] of the degen- 
erate Calogero-Moser systems and the quadratic algebras f° r the quantum systems with 
the super-integrable rational (l/(? 2 ) potential. The present article deals with one of such 
properties of the degenerate Calogero-Moser systems and it is in fact very closely related to 
the algebraic linearisation ^H]- Secondly, the nature of the multi-particle interactions of the 
Calogero-Moser systems is recognised to be governed by the root systems associated with 
finite reflection (Coxeter/Weyl) groups [TJl ESI- The original models are all based 

on the A type root system related to the symmetric group Sn with N being the number 
of the particles. The Sn is also the Weyl group of the special unitary group SU (N) . The 
integrability (the Lax pair) of the systems based on the classical root systems (A, B, C and 
D) is noticed immediately by Olshanetsky and Perelomov ^21 G3] but the actual demon- 
stration of the integrability of the Calogero-Moser systems based on the exceptional [5J E] 
and non-crystallographic root systems || took longer years. The classical universal Lax pair 
applicable for all types of potentials and for any root system [S] and the quantum universal 
Lax pair applicable for all degenerate potentials and for any root system ^3] have been 
known for some years. 

Let us denote by A a root system of rank r. It is a finite set of R r vectors which is 
invariant under reflections in the hyperplane perpendicular to each vector in A. A reflection 
s p in terms of a root p is defined by 

s p (x) = x — p (p v ■ x), x e R r , (2.1) 

in which p v = 2p/p 2 . Thus A is characterised by 

s p {r))eA, \/ P , v eA. (2.2) 

The dynamical variables are the coordinates qi G R, i = l,...,r and their canonically 
conjugate momenta Pi 6 R, i — 1, ...,r, except for the A r case in the ordinary embedding, 
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in which the number of particles is r + 1. The Hamiltonian for the classical Calogero-Moser 
system with the rational potential but without the harmonic confining potential is: 

in which the real and positive coupling constants g\ p \ are defined on orbits of the corre- 
sponding Coxeter group. That is, for the simple Lie algebra cases g\ p \ = g for all roots in 
simply-laced models and g\ p \ = gL for long roots and g\ p \ = gs for short roots in non-simply 
laced models. In order to define Lax pair matrices L and M, let us choose a set of R r vectors 
H = {a, /?,...,}, #72 = V, permuting under the action of the reflection group: 

s p (a) eTZ, VaG TZ, Vp e A. (2.4) 

We demand that it consists of a single orbit of the Coxeter group, for irreducibility. Then 
we define DxD matrices indexed by the elements of 1Z: 

p ■ H : (p ■ H) aP = (p • a) 5 al3 , (2.5) 
s P : {s p ) a p = 6 a>SpW . (2.6) 

Introduce next the T> x T> matrices X, L and M [3IH1 : 

X = i g g^p.H)-J—^ Sp , (2.7) 

L = p-H + X, (2.8) 



2 



and a diagonal matrix: 



Q = q-H: {Q) a p = {q-a)5 aP . (2.10) 

Here L and Q are hermitian = L, = Q and M is anti-hermitian = —M. 

As shown in ^U] (12. 7a), (12. 7b) the time evolution of the matrix L along the flow of the 
Hamiltonian (|2.3|) displays the following equations: 

dL . , 

-ft = i L ' M l (2-11) 

^ = [Q,M]+L. (2.12) 
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Next let us define another T> x D unitary matrix U(t) by the linear equation and the 
initial condition: 

^ = UM, U(0) = l v , (2.13) 
in which 1© is the DxD unit matrix. The final step is the introduction of W(t): 

W(t) = U{t)Q(t)U-\t), (2.14) 

which has a simple time evolution 

W = U(Q -[Q.M^U- 1 = ULU~\ (2.15) 

W = U{L-[L,M])IJ- 1 = 0. (2.16) 

The solution is 

W(t) = W(0)+tW(0), (2.17) 

with the initial values 

W(0)=Q{0), W(0) = L(0), (2.18) 

which are determined by the initial values of the canonical variables qj{0), Pj{0), j = 1, • • • , t. 

Due to the defining relation of W(t) in terms of the diagonal matrix Q(t) ()2.14j) . the solution 

{q{t)} of the the canonical equations of motion 

% = &H dpj = &H ■ = 1 r , 2 19) 

dt dpj ' dt dqj ' 

with the above Hamiltonian ()2.3|) . is simply obtained by diagonalising the above matrix 

solution ()2.17|) . The conjugate momenta {pit)} are obtained by differentiation Pj(t) = 

dqj(t)/dt. 

As promised, this is the universal proof applicable for any root system including the non- 
crystallographic one. The spectrum of W(t) ()2.17|) is highly constrained, since its dimension 
T> is usually much greater than the degree of freedom r. The high symmetry of the spectrum 
is guaranteed by the Coxeter invariance of the theory: 

H(s p (p),s p (q)) = H(p,q), VpGA, (2.20) 
L(s p (p),s p {q)) = s p L(p,q)s p , M(s p (q)) = s p M{q)s p . (2.21) 

The original proof of the explicit integration of the A type systems by Olshanetsky and 
Perelomov jH] is the very special case in which the spectrum of W(t) ()2.17|) is not constrained. 
Our proof reduces to that of jH] when A = A r and the set of vector weights is chosen as 
n = V, #K = V = r + 1. 
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3 Rational with Harmonic Confining Potential 

The Hamiltonian is now: 

With the same matrices introduced above in the preceding section, the time evolution dis- 
plays (13. 2a), (13. 2b): 

L = [L,M]-uj 2 Q, (3.2) 
Q = [Q,M]+L. (3.3) 

With the same definition of the unitary matrix U(t) as above (|2.13|) . the matrix 

W{t) = U{t)Q{t)U-\t), (3.4) 

evolves harmonically in time: 

W = U(Q -[Q.M])^ 1 = ULU- 1 (3.5) 

W = U(L — [L, M])f/ _1 = —u 2 W. (3.6) 

The solution is 

W{t) = cos cut W (0) +LJ- 1 smut W{0), (3.7) 

with the initial values 

W(0)=Q{0), W(0) = L(0). (3.8) 

Again the explicit solution {q(t}} is obtained by diagonalising the above matrix W{t) ()3.7j) 
with the harmonic time dependence. 

4 Trigonometric Potential 

The Hamiltonian of the trigonometric (Sutherland) model [2j writes: 

2 2 P tA + sm (P ■ <?) 

In order to get the hyperbolic case it suffices to change sin into sinh. In the following, we 
only demonstrate the explicit integration of the trigonometric case. The hyperbolic case can 
be deduced easily by the above replacement. 
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Two types of Lax pairs are known j^l IE] for the trigonometric cases: the minimal and the 
universal Lax pairs. While the latter, the universal lax pair, applies to any crystallographic 
root system, the former, the minimal Lax pair, requires 1Z to be the set of minimal weights, 
which exists only for the A, D, E e and E7 root systems. Let us start with the minimal Lax 
pair which has simpler structure thanks to the restriction to the minimal weights, satisfying 
the condition: 

fi : minimal weight <=> a y ■ fi = 0, ±1, Va G A. (4-2) 

4.1 Minimal Lax pair 

We consider the matrices |UJ IE] : 

L = p-H + X, (4.3) 
X = iJ^(,.^_L_*,, (4.4) 

m = -\ e J p{2c :; {p -f is p -i v)+ i e ^-fe^c, (4.5) 

2 P tX + sm (p ■ 9) pl^ + sm (p • 9) 

and diagonal matrices: 

Q = q-H : (Q) Qj g = (g • a) (4.6) 

R = e 2iQ . (4.7) 

Again L and Q are hermitian = L, = Q and M is anti-hermitian M' = —M. Thus i? 
is unitary. 

As shown in (TUj (15. 3a), (15. 3b), when the root system admits a minimal representation, 
and 1Z being the set of minimal weights, the time evolution along the flow of the Hamiltonian 
(I4~H displays: 

r)T 

-ft = [L,M], (4.8) 
ri n 

— = [R, M]+i (RL + LR) . (4.9) 
With the same definition of the unitary matrix U(t) as above (|2.13|) . 

^ = UM, U(0) = l v , (4.10) 
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we introduce a matrix 

W(t) = U(t)R(t) Uit)- 1 = U(t) e 2iQ{t) U(t)-\ (4.11) 

It satisfies a simple first order linear differential equation 
dW 



Of 



U(dR/dt - [R, M})^ 1 = i U(RL + LR)^ 1 (4.12) 
ifVVULU- 1 + ULU~ 1 W) , (4.13) 



since as in (|2.15|) . (|2.16|) . ULU 1 is a constant matrix: 

d_ 

Ft 



a 

— (UL U~ l ) = U{dL/dt - [L, M])f/ _1 = 0, (4.14) 



UtyLtyUit)- 1 = L(0). (4.15) 

The solution is 

By diagonalising the above matrix solution, we obtain the explicit solution {q(t)} of the 
classical Sutherland model (|4.1|) . One might naturally wonder if the coordinates {q(i}} could 
be determined uniquely from the unitary matrix (|4.16jl . The answer is affirmative since the 
motion is always restricted to one of the Weyl alcoves due to the periodicity and singularity 
of the potential. Near the boundary of a Weyl alcove, for example at p ■ q = 0, p G A, the 
singularity of the potential ~ l/(p-q) 2 can never be surpassed classically. Therefore if {q(0)} 
is in the principal Weyl alcove, 

PW T = {qER r \p-q>0, pEll, p h -q<n}, (4.17) 

{q{t)} will always remain there. Here II is the set of the simple roots and ph is the highest 
weight. This removes any ambiguity in determining {q(t)} from the eigenvalues of ()4.16|) . As 
in the rational potential cases, the spectrum of W(t) is highly constrained as a consequence 
of the Weyl invariance IjOljl . (l2~2H . 

4.2 Universal Lax pair 

The universal Lax pair has cot(p • q) function in L instead of 1/ sin(p • q) in (14.41) 

L = p-H + X, (4.18) 

X = i E 9\ P \(p-H)cot(p-q)s p , (4.19) 

M = -- E 9\ P \ ■ if , (s P ~ Id) , (4.20) 
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which satisfy dL/dt = [L,M] for the Hamiltonian flow but the additional equation ()4.9|) 
takes a different form. 

For 1Z being the set of minimal weights, it reads 

— = [R J M}+i (R(L + K) + (L- K)R) , (4.21) 
in which K is a non-negative constant matrix commuting with M: 

K= £ 9\ P \(P-H)(p v -H)s p , [K,M} = 0. (4.22) 

It is a very important quantity in Calogero-Moser systems appearing in many contexts. For 
example, it is a commutator of Q (|2.1(Jj) and the rational Lax matrix L (|2.8|) . (J2.7|) (see (4.36) 
of US] and (2.40) of [E]): 

[Q, L] = %K. (4.23) 

It should be noted that if K is defined as above, the expression (J4.22|) is universal, that is 
valid for any root system A and any choice of TZ. Various properties of the K matrix, whose 
eigenvalues are all 'integers', are discussed in detail by Corrigan-Sasaki, in the Appendix of 

For TZ being the set of all roots A (for the simply-laced root systems) or the set of short 
roots As (for non simply-laced root systems) and also the set of vector weights (V) for the 
C, the relation corresponding to (|4.9|) and (|4.21|) reads 

— = [R, M] + i (R(L + K) + (L- K)r) , (4.24) 
in which K is another constant matrix commuting with M, 

k = E 9\ P \\P-H\s p , [K,M}=0, (4.25) 

introduced by Corrigan-Sasaki, as (5.32) of [To] . 

Now the explicit solution of the Sutherland system is achieved for any crystallographic 
root system, since one can choose at least one Lax pair satisfying (|4.24j) . We proceed as 
before by defining the unitary matrix U(t) by (j4.1(Jj) and introduce a matrix 

W(t) = U(t)R(t) Uit)- 1 = U(t) e 2lQit) U(t)-\ (4.26) 
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It satisfies a simple first order linear differential equation 



U(dR/dt - [R, M])U~ l = i U{R{L + K) + (L - K)R)U 

i (w u(l + k) tr 1 + u(l - k) tr 1 w) , 



-i 



(4.27) 



dt 



(4.28) 



since as in (j2~T5j) . (|2~TT)J) . U(L ± K) U' 1 i 



is a constant matrix: 



-(U(L±K)U- 1 ) = U(dL/dt-[L,M])U- 1 = 0, 
U(t)(L(t) ±K)U(t)- 1 = L(0)±K. 



(4.29) 
(4.30) 



The solution is 



Wit) = e 



,it(L(0)-K) e 2iQ(0) e it(L(0)+K) 



(4.31) 



By diagonalising the above matrix solution, we obtain the explicit solution {q(t)} of the 
classical Sutherland system (|4.1j) for any root system. 

The very fact that K (K) commutes with M simply means that a spectral parameter A 
can be introduced trivially into the Lax pair for degenerate potentials jH] : 



5 Rational Higher Flows 

The integrable hierarchy of the Calogero- Sutherland systems consists of Hamiltonians gen- 
erated by higher conserved quantities, which are constructed, for example, from the trace 
of the higher powers of the L matrix, 7i n oc Tr(L 2n ). The method of explicit integration as 
described in the preceding sections applies also to these higher Hamiltonian flows, as shown 
by Suris [7] for the A type root systems with the conventional Lax pair, 1Z = V. However, 
in contrast to the basic Calogero-Sutherland flows, it works only for those systems based on 
the classical root systems, the A, B, C and D. Let us denote by 1Z the set of vector weights 
V, for the A, C and D root systems and the set of short roots As for the B root system. 
These particular sets 1Z have unique orthogonality property 



which endows a very special structure to the Lax pair represented on 1Z. The dimensions of 
the corresponding Lax matrices are T> = r + 1 for the A r and T> = 2r for the B r , C r and D r . 



L X = L + \K, L X = [L X ,M}. 



(4.32) 



if fi 7^ ±u, fi ■ v = 0, V/i, v G 7Z, 



(5.1) 
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It is through these special Lax matrices that the explicit integration of the higher rational 
and trigonometric flows is realised. 

Let us start with the explicit forms of the rational L matrices: 

(A) : L, L jk = pj5 jk + ig(l - 5 jk )/(qj - q k ), (5.2) 

(B) : L =( K _ B ' A jk =p j 5 jk + ig L (l- 5 jk )/(qj - q k ), 

B jk = i(g s /qj)5 jk + ig L (l - $jk)/(qj + qk)- (5.3) 

The rational C system will not be discussed since it is equivalent to the rational B system. 
The rational D system is obtained by constraining g$ = in the rational B system. 
The higher Hamiltonians are 

(A): H n = Tr(L n+1 )/(n + l), n > 1, (5.4) 
(B,D): H n = Tr(L 2n )/(4n), n > 1. (5.5) 



The lowest 7i\ is the original Hamiltonian ()2.3)1 . The basic idea is to rewrite the Hamiltonian 
flow 

dqj &H n dpj dH n 



dt n dpj ' dt n dqj 

into equivalent matrix forms 



(5.6) 



df n = [L ' Mn] ' (5 - 7) 
|g = [Q,M n ]+L n (L 2 - 1 ), (5.8) 

as in the lowest flow (l2~TTJ) . (ETHl . 

In contrast to the lowest flow case in which the explicit form of M is given ()2.9|) . we can 
interpret part of ()5.7|) and ()5.8|) as determining M n . The diagonal part of the Q equation 
()5.8|) is equivalent to the first half of the canonical equations (15. 6|) . The off-diagonal part of 
the Q equation ()5.8|) determines the off-diagonal part of M n completely: 

(M n ), u = -(L"W<? -(ji-v) (-(L 2 "- 1 Wg •(//-»/)), \i±v. (5.9) 

Whereas the diagonal part of M n does not enter the the Q equation ([5.8)1 . it can be deter- 
mined from the off-diagonal part of the Lax equation. The result is very simple: 

(M n ) w = - J2 {M n )»n = - E ( M «W M l = -M n . (5.10) 
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The proof that the diagonal part of the higher flow Lax equation ()5.7|) is equivalent to the 
second half of the canonical equations ()5.6|) goes almost parallel to that of the lowest flow. 

After the equivalence of the canonical equations (|5.6|) with the two matrix equations ()5.7|) 
and ()5.8|) is established, the explicit integration by diagonalisation is straightforward. Let 
us define aDxD unitary matrix U n (t n ) by the linear equation and the initial condition: 

8U 

—^ = U n M n , U n (0) = l v . (5.11) 
ut n 

Then a matrix function W n (t n ), defined by 

W n (tn) EE U n (t n )Q(tn)U^(t n ), (5.12) 

has a simple time evolution 
dW n 



dt n 
_d 
dt, 

The solution is 



U n (dQ/dt n - [Q, M n \)U- 1 = {U n LU- x ) n ((UnLU- 1 ) 211 - 1 ) , (5.13) 

° (u n LU- 1 ) = U n (dL/dt n - [L,M n ])U? = 0, (5.14) 

U n (t n )L(t n )U n (t n )- 1 = L(0). (5.15) 

W n {t n ) = W n {0) + t n dW n {0)/dt n , (5.16) 



with the initial values 

W n (0) = Q(0), dW n (0)/dt n = L(0) n (L(O) 2 "- 1 ), (5.17) 

which are determined by the initial values of the canonical variables qj{0), Pj{0), j = 1, • • • , r. 
Due to the defining relation of W n (t n ) in terms of the diagonal matrix Q(t n ) (j5.12)l . the solu- 
tion {q(t n )} of the the canonical equations of motion ()5.6j) with the above Hamiltonian ()5.4j) 
or ()5.5jl . is simply obtained by diagonalising the above matrix solution (|5.16jl . Determina- 
tion of the conjugate momenta {p{t n )} requires solution of the second half of the canonical 
equations of motion (|5.6J1 . which are now algebraic since {dq/dt n } are now known functions 
of time. Extension to the generic higher flows of the hierarchy 

H = '^2c n TC n , c n : const., (5.18) 

n 

is straightforward since the matrix equations ()5.7|) and (|5.8J) are linear in M n . However, 
some higher flows cannot be treated this way. For example, in the D r (r: odd) theory, there 
exists another conserved quantity (Hamiltonian) of the form pip 2 ■ ■ ■ p r + ■ ■ ■, which cannot 
be written as (|5.18|) . 
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6 Trigonometric Higher Flows 



The basic logics of the explicit integration of the trigonometric higher flows is almost the 
same as that of the rational higher flows, except that we have to consider two different types 
of Lax pairs; the minimal and the universal. So we just write down the key formulas without 
detailed derivation. 

6.1 Minimal Lax Pair 

We discuss the explicit integration of the trigonometric higher flows of the A and D theory 
in terms of the minimal Lax pair, although the formulation in terms of the universal Lax 
pair works well for them, too. 

The explicit forms of the trigonometric minimal L matrices are: 

(A): L, Lj k — pj5j k + ig(l — 5 jk )/ sin(^- — q k ), (6.1) 

(D) : L=(^_ B _ A \ A jk =p j 5 jk + ig(l-5 jk )/sm(q j -q k ), 

B jk = ig(l - 8 jk )/ sin(g j + q k ). (6.2) 

The higher Hamiltonians take exactly the same form as ()5.4j) and (|5.5j) . The lowest 7i\ is 
the original Hamiltonian (|4.1j) . We rewrite the higher Hamiltonian flow ()5.6j) into equivalent 
matrix forms 

dL 

Of = [L,M n ], (6.3) 

a d 

= [R, M n ] + i(RL n + L n R) (i(RL 2n ~ l + L 2n ~ l R)^j , (6.4) 
as in the lowest flow ()4.8|) . ()4.9|) . The off-diagonal part of M n is 

(M n )^ = -(L n )^ cot[g ■ (ji - u)] (-(L 2 "" 1 )^ cot[g ■ (/i - u)]) , fi ^ v. (6.5) 
The diagonal part is 

(M n ) w = (L n U/sm[q-(iy-v)} = £ {L n )^ v / sin[q • (// - v)], = -M n . (6.6) 

The VxD matrix W n (t n ) obeys simple time evolution: 

W n (t n ) = U n {t)R(t n )U n {t n )- 1 =U n (t n )e 2iQ ^U n (t n )-\ (6.7) 

e it n L(0) n e 2iQ(0) e it n L(0) n ^„L(0) 2 "- 1 e 2iQ(0) e it„L(0) 2 ™- 1 ^ _ 
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By diagonalising the above matrix solution ([6.8)1 . we obtain the explicit solution {q(t n )} of 
the higher flows of the Sutherland system (|5.4|) and ()5.5|) for the A and .D root systems. 

6.2 Universal Lax Pair 

The explicit integration of the higher flows of the B and C Sutherland systems is achieved 
in terms of the universal Lax pairs based on the set of short roots (1Z = As) for B and the 
set of vector weights (TZ = V) for C. For the rank r system both have T> = 2r. 
The Lax matrix L and the constant matrix K (j4.25j) are: 

l =(-b -a)' Mrs)' (6J ) 

(B) : A jk = pj8 jk + ig L {l - S jk ) cotfy - q k ), S jk = g L (l - 5 jk ), (6.10) 

B jk = ig s cot qj 5 jk + ig L (l - 5 jk )cot(qj + q k ), T jk = g s 8 jk + g L (l - S jk ), (6.11) 

(C) : A jk = pj5 jk + ig s (l - 5 jk ) cot(g j - q k ), S jk = g s (l - 5 jk ), (6.12) 

Bjk = 2ig L cot2q j 5 jk + ig s (l - 5 jk ) cot(g j + q k ), T jk = 2g L 6 jk + g s (l - 5j k ).(6A3) 

It is easy to see 

e iQ {L + K)e~ iQ = e~ iQ {L - K)e iQ , (6.14) 
^Tr(L + K) n = Tr(L — K) n , (6.15) 

which are conserved quantities of the Sutherland flow (|4.ip . It differs from the usual one 
Tr(L n ) by a linear combination of lower order conserved quantities. The canonical equations 
of the higher flow Hamiltonian 

H n = Tr((L±K) 2n )/(4n) (6.16) 

are equivalent to the matrix equations 
dL 

df n = [L ' Mn] ' (6 - 17) 

— = [R,M n }+i(R(L + kf n - l + (L-k) 2n ~ l R). (6.18) 

The off-diagonal part of M n is 

(M n ) MV = - [e^-»\L + Kf^ + e -^- v \L - K) n J /sin[g ■ (ji - u)}, fi ^ v. (6.19) 
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The diagonal part is 



(M n ) 



M\ = -M n . 



(6.20) 



The D xP matrix W n (t ri ) obeys simple time evolution: 



W n (t n ) 



U n (t)R(t n ) U n (t n ) 

e ii„(L(0)-K) 2 "- 1 e 2i(; 




(6.21) 
(6.22) 



By diagonalising the above matrix solution ()6.22|) . we obtain the explicit solution {q(t n )} of 
the higher flows of the Sutherland system ()6.16j) for the B and C root systems. 

7 Summary and Comments 

Explicit integration of the Calogero and Sutherland systems by means of diagonalisation 
is demonstrated for any root system, the exceptional as well as the classical and the non- 
crystallographic. It is based on the universal Lax pair for the degenerate potentials, that is 
the rational with/without the harmonic confining potential and the trigonometric/hyperbolic 
potentials. As emphasised in the text, it is very closely related to the concept of algebraic 
linearisation by Caseiro-Frangoise-Sasaki ^0] • The method is extended to the higher Hamil- 
tonian flows of the rational and trigonometric/hyperbolic interactions. In contrast to the 
basic Calogero-Sutherland flows, the applicability is limited to those systems based on the 
classical root systems, the A, B, C and D root systems. 

The theory of explicit integration of higher Hamiltonian flows is very closely related to 
the dynamical r- matrix and the Hamiltonian reduction fOlEj- In the case of the 

most classical rational potential of the A type, the method of Hamiltonian reduction starts 
from the large phase space of the matrix dynamical variable W and its conjugate momentum 
variable Z, which are both assumed to be Hermitian. The Hamiltonians 



H n = Ti(Z n+1 )/(n + l) 



(7.1) 



generate the flows 



dW 



dZ 



(7.2) 



0. 
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This Hamiltonian system in invariant under the action (W,Z) — > (UWU 1 ,UZU l ) of 
unitary matrices U. The reduced phase space is obtained by imposing the constraint 



[W, Z] = iK 



(7.3) 



and factoring out the constrained phase space by residual symmetries (i.e., by the group of 
unitary matrices that commute with K). The (Q, L) pair ()4.23j) is nothing but a representa- 
tive of a point of the reduced phase space, which is connected with the point (W, Z) of the 
large phase space by a (t-dependent) unitary matrix U as 



The linear flows of (W, Z) are thereby mapped to the Calogero flows of (Q, L). This is the 
way to understand the rational Calogero system of the A type as a Hamiltonian reduction 
19 ; a similar interpretation has been proposed for a few other cases [T3]. The dynamical 
r-matrix has been constructed in this framework of Hamiltonian reduction jjj. We expect 
that all the cases discussed in this paper can be treated in the same way. 
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